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SUMMARY 


A linear system for applying thrust to a ferry vehicle in the 
terminal phase of rendezvous with a satellite is analyzed. This system 
requires that the ferry thrust vector per unit mass be variable and 
equal to a suitable linear combination of the measured position and 
velocity vectors of the ferry relative to the satellite. The variations 
of the ferry position, speed, acceleration, and mass ratio are examined 
for several combinations of the initial conditions and two basic control 
parameters analogous to the undamped natural frequency and the fraction 
of critical damping. Upon making a desirable selection of one control 
parameter and requiring minimum fuel expenditure for given terminal- 
phase initial conditions, a simplified analysis in one dimension prac- 
tically fixes the choice of the remaining control parameter. The system 
can be implemented by an automatic controller or by a pilot. 


INTRODUCTION 


The problem of the rendezvous of space vehicles is one of increasing 
significance in space research and engineering. References 1 to 14, 
which are typical of the available literature in this field, contain 
analyses of various aspects of the rendezvous problem. In these papers 
attention has been given mainly to the problem of most immediate interest, 
namely, that of the rendezvous of an earth- launched ferry vehicle with a 
satellite or space station in orbit about the earth. 

From the standpoint of thrust application, terminal- phase rendezvous 
systems include (a) those which utilize one or more constant- thrust rocket 
motors, and (b) systems with continuously variable thrust. Terminal- phase 
rendezvous systems employing continuously variable thrust have been ana- 
lyzed. (For example, see refs. 8, 13, and 14. ) Systems of type (b), 
another of which is considered in this paper, can be approached in prac- 
tice by systems of type (a); thus, the throttlability limitations of 
existing rocket motors may not be a serious obstacle to the design of a 
system which, in its simplest form, would utilize one or more throttlable 
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motors. Methods for approaching a system of type (b) by systems of 
type (a) are illustrated schematically in figure 1 . 

In this paper, thrust which is continuously variable in magnitude 
and direction will be assumed to be available. The basic principle of 
the thrust control system to be employed requires that the instantaneous 
vector acceleration of the ferry be made equal to the sum of two vectors, 
which are suitable multiples of the position and velocity vectors, 
respectively, of the ferry relative to the satellite. Thus, the accelera- 
tion vector of the ferry is made to be a suitable linear combination of 
its relative position and velocity vectors. The three vectors are then 
coplanar. For the analysis of this system, a nonrotating satellite- 
centered coordinate system will be used. (See fig. 2. ) The choice of 
the control parameters, that is, determining the linear combination of 
the position and velocity vectors that shoula be used, is examined herein. 
Moreover, the effects of the initial conditions and of gravitation on the 
terminal-phase motion and the rocket fuel expenditure are discussed. 


SYMBOLS 


a 


F 

G 


m 

mo 

M 

R 

t 

T 

Vi,V 2 


distance between satellite and ferry at start of terminal 
phase of rendezvous 

ferry thrust vector 

Newton's universal gravitational constant, 

6 . 67 O X 10“^ newton-meter^/kg^ cr 3.458 X 10"® lb-ft^/slug^ 

orthogonal unit vectors 
ferry mass 

ferry mass at start of terminal phase of rendezvous 
limit of m as t approaches 00 

mass of earth, 5*975 X 10^ kg or * .094 X 10^ slugs 

distance between satellite and feny 

time 

satellite orbital period 

ferry relative velocity components at start of terminal 
phase of rendezvous 


h ir\ tc\ lov 
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r l> v 2 


v ( 


ex 


nondimens ional values of ( re ^ erre( i to u\) a ) 

effective exit speed of propulsive exhaust gases 

ferry position vector in a coordinate system with origin at 
satellite and axes always parallel to lines fixed in an 
inertial frame 


x l’ x 2’ x 3 com P° nen ^ s vector x 


(X~i 


value of x-l at t -- 0 


X 


control parameter analogous to fraction of critical damping 

ferry position vector in a coordinate system with origin at 
center of earth and axes oriented as for x above 


-1 x ? 

angle , tan -=• 
X 1 


mass variable, 


m 


v ex/^O a 




£ satellite position vector in a coordinate system with origin 

at center of earth and axes oriented as for x* 

£ distance between satellite and center of earth 

^1* ^2*^5 com P° nen t s of vector X 


O3 0 


control parameter analogous to undamped natural frequency 


Cb 


= when £ < 1 


Q = o) 0 ^ 2 - 1 when £ > 1 


Dots over symbols indicate differentiation with respect to time. 
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SOME GENERAL CONSIDERATIONS ON THE TERMINAL PHASE OF RENDEZVOUS 


In order to have a basis for ideas, it nay be useful to define the 
start and end of the terminal phase of rendezvous as the times at which, 
respectively: 

(a) Radar or other contact has been established between ferry and 
satellite to permit measurements of relative position and velocity and 
based on these measurements, final- approach Ihrusting is initiated to 

reduce the relative distance and velocity to low values (for example, p 

100 feet and 5 feet per second). q 

(b) Relative position and velocity have been reduced to sufficiently 
low values to permit the initiation of docking. 

Before proceeding to the discussion of the thrust-application system 
which is the subject of this paper, brief consideration will be given to 
some simple concepts which are useful in rendezvous studies when the 
effect of gravitation on the motion of the ferry relative to the satel- 
lite is a minor one. 

If the velocity of one vehicle is measured with respect to a stable 
platform in the other vehicle, the analysis presented in reference 15 
(pages 5 to 10) can be applied directly to the motion of the ferry rela- 
tive to the satellite provided the gravitational field is approximately 
central. In particular, the analysis indicates that gravitation may be 
of secondary importance (depending upon the thrust level, the apparent 
gravitational acceleration, and so forth) in its effect on the relative 
motion in the terminal phase of a given rendezvous situation. If this 
is the case, rendezvous systems can be subjected to approximate analyses 
which neglect the gravitational effect. The limitations of simplified 
studies of this type can subsequently be appraised in a number of ways, 
for example, by computing the time integral of the apparent gravitational 
acceleration and comparing this integral with, a speed characteristic of 
the ferry's approach. In this way, more detailed and accurate studies 
of trajectories, fuel consumption, and so foith may be deferred until 
the final stage of a design program. 

Consider a nonrotating coordinate system (see fig. 2) with origin 
at the satellite (target) and assume that gravitational effects on the 
relative motion are negligible. In this reference frame, the primary 
forces acting on the ferry are the thrust anc orientation-control forces. 

The ferry speed (relative to this reference frame) at the start of the 
terminal phase puts a lower bound on the required fuel expenditure. 

Rendezvous with this minimum expenditure of fuel can clearly be accom- 
plished in a finite time only if the ferry relative velocity vector 
points directly toward the origin (the posit:’ on of the satellite in the 
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reference frame). In general, then, it is advantageous to start the 
terminal phase of rendezvous with the ferry velocity vector (l) directed 
as nearly as possible toward the origin, and (2) of just sufficient mag- 
nitude (speed) to permit rendezvous to be accomplished in a specified 
time without the necessity of extra expenditure of fuel to hasten the 
closure . 

For example, in a nonrotating, earth-centered reference frame, a 
ferry may be launched into an elliptical orbit which osculates the 
satellite orbit from inside the latter. (See fig. 3(a)- ) It the flight 
is planned so that the ferry arrives at the oscule (osculation point) 
somewhat ahead of the satellite and so that the terminal- phase thrusting 
begins at this time, the expenditure of fuel in the terminal phase prin- 
cipally completes the process of bringing the speed of the ferry up to 
the speed of the satellite. On the other hand, in the case of a ferry 
orbit which osculates the satellite orbit from outside the latter 
(fig. 3 (b) ), it would be preferable to have the satellite reach the 
oscule ahead of the ferry, at which time the terminal- phase thrusting 
is initiated and is utilized mainly to bring the speed of the ferry down 
to the speed of the satellite. In both examples, small variations of 
the direction of the thrust vector which might be required for rendezvous 
maneuvering would add little to the fuel expenditure. However, the sys- 
tem to be analyzed is not restricted to these situations. 


A LINEAR SYSTEM FOR THRUST APPLICATION 


The basic requirement for the terminal approach, the near nullifi- 
cation of the relative position and velocity vectors, calls for measure- 
ments of either these vectors or an equivalent set of variables. 

One of the simplest types of terminal- phase rendezvous systems 
utilizing continuously variable thrust is one in which the thrust vector 
per unit mass is made a linear combination of the relative position and 
velocity vectors of the ferry. (See fig. 2. ) The following analysis 
of this type of system at first neglects the effect of gravitation on 
the relative motion. Subsequently, the central-field gravitational 
effect is examined in some cases of current interest. Trajectory per- 
turbations due to drag, electric or magnetic fields, solar radiation 
pressure, and so forth are neglected throughout the paper. Similarly, 
problems concerned with the orientation of the ferry as a rigid body 
are not considered. 

In a rectangular Cartesian coordinate system with origin at the 
satellite and with axes which maintain constant directions with respect 
to an inertial frame, that is, with respect to the fixed stars, let the 
position vector of the ferry be denoted by 
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x li + 


x 2 J 


x^k 


(i) 


where i, j, and k are unit vectors along the coordinate axes. It 
should be noted that the negative of x is tae position vector of the 
satellite with respect to a stable platform or a star-oriented system 
of axes in the ferry. Under the simplifying assumptions of the preceding 
paragraph, the vector equation of relative motion is 


, 2 — » 

d x _ F 


dt 


m 


( 2 ) 


where m is the mass of the ferry and the vector F is the thrust 
acting on it. 


In the terminal-phase rendezvous system considered here, thrust is 
applied in accordance with the equation 


F 

m 




(3) 


where cdq and £ are suitably chosen constants. If the right-hand 

side of equation (j) is considered as the desired thrust per unit mass 
and the left-hand side as that which is actually applied, it is clear 
that the equation can be only approximately satisfied in practice. The 
action which equation ( 3 ) calls for can be effected in a number of ways 
by an automatic controller or by a pilot. The limiting case of a perfect 
controller is considered in the analysis that follows, that is, equa- 
tion ( 3 ) is assumed to be exact. 


Equations (2) and ( 3 ) yield 


d x 
dt 2 


2 ^0 if + 


-x = 0 


(M 


which is the equation of motion of the ferry "elative to the satellite 
for the system under consideration. Equation (4) is equivalent to the 
equation of motion of a harmonic oscillator in three dimensions with 
damping proportional to the velocity. The second term, or damping term, 
in equation (4) may be thought of as the "slowing" part, whereas the 
third term, or restoring term, could be called the "zeroing" part. The 
latter, however, is only a tendency. Except in the trivial case with 
the initial conditions (at t = 0) 


x = 0 


k] H LOi ro IA 
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equation (4) yields no solution corresponding to rendezvous in the 
precise sense 

5 ? = 0 



(at t = Some finite time) for these final conditions imply x - 0 for 
all t (including t = 0); which is the trivial case mentioned. The 
proof; which is a simple consequence of the theory of linear differential 
equations with constant coefficients; is omitted. 

Rendezvous in the practical sense; that is; the near nullification 
of x* and dx/dt within a reasonable time; can be accomplished with 
the system described by equation (4). Moreover; the system possesses 
certain desirable features besides the simplicity of its analytical 
expression. These; together with the most questionable feature; the 
assumption of continuously variable thrust; will be discussed subsequently. 

It follows from its vector nature that equation (4) is invariant 
under rotations of the coordinate system. For convenience; then; let 
the unit vector i have the direction and sense which the position 
vector x has at t = 0; the start of the terminal phase. Thus; at 
t = 0 


x=ai (a>0) (5) 

It is also convenient to choose the unit vector j so that (see fig. b) } 
at t = 0; 

= V x i + V 2 T (v 2 £ o) (6) 

This choice of j* is evidently unique unless V 2 = 0. In this case 

the motion is confined to the line of i*. Otherwise; the trajectory 
determined by equations (4); (5); and (6) is confined to the plane of 
1 * and j; that is, x^ = 0 for all t. 

The parameters cdq and £ in equation (k) are easily recognized 
as the undamped natural frequency and the fraction of critical damping; 
respectively; of the problem of the damped harmonic oscillator. Thus; 
multiplication by ooq renders t dimensionless. 



Solution for £ < 1 


The solution of equation (4) subject to :he initial conditions ( 5 ) 
and (6) can, for £ < 1, be expressed in the dimensionless form 


i - [cos ot * (5 ♦ v, )2±2^y 

8 1 ^ifiTT 2 


x 2 sin art ~ 

-5- = v 2 y == e 


W - 


1 ~ a) 0 a 


V 2 o^a 


The corresponding velocity components are giv?n by 


= fv. cos art - (l + 

“o a 1 1 1 \r~"72\ 


yi - 


x 2 / , „ sin cut \ -^owt 

= ( v 0 cos cot - r v 0 ] e u 

“° a \ \f - 


where the dots refer to differentiation with aspect to time. The 
acceleration components are, then, 


«\> 2a 


-(1 + 2 5 v,)co S a* + [? + (2c 2 - 

V 1/ N ' - 1 I / O 1 


2 

cjoq a 


= -2^v 2 cos cut + (2£ 2 - l)' 


sin cot ; - £co n t 


vji \j] h tr 1 
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Equations (7), ( 9 ), and (10) correspond to the underdamped (£ < l) 
harmon i c o s c i llato r . 


L 

L 

3 

3 

3 


Solution for £ = 1 


For critical damping (£ = l), the characteristic equation of equa- 
tion (4) has equal roots. In this case, the solution of equation (4), 
again subject to equations (5) and (6), takes the form 


-1 


a 
x 2 

— = v 2 o) 0 te 


1 + (1 + vi)o) 0 tj e 0 


"'N 


-CUQt 


( 11 ) 


with velocity components 


tDQa 


Vq_ - (l + v 


"\ 
-(JD n t 


fV] 6 " 0 


-Mrs t 


= v P ( 1 - owt)e 0 
o)Qa ^ \ U / 


(12) 


and acceleration components 
x 


k l 


p 

^0 a 


2 

o)Q^a 




-(! + 2v x ) + (1 + v-^tje ° 


= vp(-2 + u>Qtje 






(13) 


Solution for £ > 1 

Finally, for the overdamped case (^ > l), the solution has the form 


X 1 , /„ > sinh Qt 

— - [cosh nt + ^ + v i) 




-^ 0 t 


?\ 


(14) 


x 2 sinh fit 

— - v Q — e u 

a ^ 


r 7 


j 
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where fi = cdq)J^ - 1. The corresponding velocity components are 


A 1 

o^a 


v^ cosh fit - (! + £ v l) 


si nl fit 


£ 


-^ 0 t 


is. . Vpfcosh at - 5 
&Q& 2 l b 


]/p 


( 15 ) 




and acceleration components are 

= f-( 1 + 2^ Vl )cosh fit + [5 + (2^ 2 - l)vi] e'^Q 1 




fi 2 - 


> (16) 


CDq' 


2 - v 2 


-2£ cosh fit + ( 25 ^ - l) — 


inh fit 


yt 2 - 1 






Mass Variation 

For any value of £ , momentum considerations yield the well-known 
rocket equation 


m 


dt = -v dm 
ex 


(17) 


where v gx is the effective exit speed of the propulsive exhaust gases. 

With the separation of variables and integration, equation (17) yields 
the dimensionless equation for the mass m of the vehicle at time t. 




(18) 


where mQ denotes the vehicle mass at t = 0. This equation will be 

used in the subsequent study of fuel consumption. The quantity on the 
left-hand side of equation (l8) is called the nass variable and it is 
independent of v gx - 


^ r| LAfAlfN 
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CRITICALLY DAMPED CASE IN ONE DIMENSION 


Before a general examination of the problem is made, a restricted 
class of cases will be considered. The earlier discussion concerned 
with figure 3 indicated the fuel- consumption advantage to be gained from 
rendezvous arranged so that the relative velocity of the ferry is directed 
as nearly as possible toward the satellite, that is, v^ < 0 and - 0 . 

Accordingly, this one-dimensional motion (in the absence of gravitation) 
will now be examined for the intermediate case where £ = 1. Fuel con- 
sumption is minimized in these cases by requiring that x’j_ not change 
sign during the terminal phase. Hence, it is of interest to focus atten- 
tion on those cases for which x-j_ ^ 0 for all t ^ 0. This will be true 
if (and only if) 

-1 S v x s - i (19) 


which follows readily from the first of equations ( 13 )- 
purpose of design, a reasonable choice is 


v, = — 


3 

k 


Thus, for the 


( 20 ) 


which is equivalent to the frequency choice 


«0 = - 


k Vq 
3 a 


(21) 


For this choice, equations (ll), (l2), and (13) reduce to 



a 


h 

Vi 


h 

a i 



( 22 ) 


where is the value of x^ at t 


0 , and equation (l 8 ) reduces to 



(23) 


Since cjdq > 0 (when Vq < 0), equation ( 23 ) becomes, in the limit as 
t — ^ 
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“0 


or, for -V^ « v ex , 


“0 


1 



where 


= lim m 
t — > 00 


(24) 


(25) 


( 26 ) 


Equations (22) are plotted in dimensionless form in figure 5* It 
is evident from the equations themselves that the proportion 

a Vl * ^1 

which initially has the value 1:1:1, approaches the value (referred to 
norm Xq/ a ) 


with increasing time; however, this ratio is cnly very roughly obtained 
within reasonable times. 


As a specific example, with £ = 1, Vq = -3/^j Vg = 0, let 
a = 100,000 feet and -Vq - 500 feet per seccnd. Then, coq = l/l50 per 
second and aq =2.2 feet per second per seccnd. Finally, at time 
t = 15 minutes, the distance has been reduced to xq = 620 feet, the 
ferry is approaching at speed -Xq= 3*7 feet q er second, and the accel- 
eration due to thrust has dropped to Xq = 0.C22 foot per second per 
second. If the rocket motor were shut down at this point, the ferry 
would coast^" to the satellite in an additional seconds or 2.8 minutes. 

j* ( 


Finally, if v ex = 10,000 feet per second, the mass loss due to fuel 

consumption is only about 5 percent of the initial mass. The magnitudes 
in this example, with the possible exception cf the 100:1 reduction in 


^The first two of equations (22) can be csed to show that the 
x 1 

coasting time — — is nearly independent of the thrust cut-off time t c - 

As t P increases from 0 to <», the coasting time decreases from - — to 
x >^0 

s* ’ 
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thrust acceleration, seem to be reasonable for the terminal phase of 
rendezvous with a near earth satellite. Moreover, the effect of the 
apparent gravitational acceleration can be shown to be of secondary 
importance in this case in both the thrusting and coasting phases. 

Since a rocket motor which can be throttled to 0.01 of maximum 
thrust is not likely to be developed within the foreseeable future, the 
implementation of the present rendezvous system for the specific case 
just discussed would require either a suitable combination of perhaps 
two to three variable- thrust motors or a scheme approaching variable - 
thrust performance. (See fig. 1.) A less extreme throttling ratio 
would result in the case considered, of course, if a higher coasting 
speed than 3*7 feet per second were tolerated. The discussion in the 
next section will clarify this statement. 

In this specific example, the system affords adequate time to pre- 
pare for docking or, in the event of some malfunction, to take emergency 
actions that may be required. Figure 5 shows that most of the relative- 
speed reduction occurs early in the thrusting phase. Therefore, since 

the relative velocity h? = x-^i* + x^j* will ordinarily not be in line 

with the origin in the more general case for which V 2 f 0, the proba- 
bility of collision at high relative speed due to rocket failure must 
be less for this system than for systems utilizing higher thrust levels 
late in the terminal phase. This property, together with some other 
matters of practical importance, will be made clearer by the more detailed 
study that follows. 


DETAILED STUDY WITH GRAVITATION NEGLECTED 


Equations ( 7 ), ( 9 ) to (l6), and (l8) have been used to compute the 
quantities 
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... / m \ Vex /“O a 

^ W 


for dimensionless times in the interval 0, l4. All combinations of the 
parameters 


▼l = -1.0, -0.5 
v 2 = 0, 0.25, 0.5 

5 =0.7, 1.0, 1.3 

were used. The values of v^ which were chosen are the extremes found 

in the one-dimensional analysis above. (See expression (19).) The 
results are shown in figure 6. The mass ratio m/mo can be found from 
the mass variable p by making use of the identity 


— = expl 
fflQ p ! 


'uy.a > 

l0 s e n 

l ex e , 


Since a)Qa « v ex for cases of practical interest, the mass ratio mytnQ 

is much nearer unity than is the corresponding value of p. For example, 
in the numerical example of the previous section, m/mo ** 0-95, whereas 

M- *» exp(-3/4) « 0.47. In the absence of gravitation, the optimal final 

(t = oo) value of p is, by equation (l8), ex;> (-f + v 2 ^ ^ • This ideal 

value can be achieved, in principle, by applying thrust to reduce the 
initial velocity to zero; infinitesimal impulses are then used to com- 
plete the rendezvous (infinite time being required unless V 2 =0). 

Consider first those cases in figure 6 for which £ = 1.0 and 
£ = I- 3- It is noted that the relative distance, speed, and accelera- 
tion decrease approximately exponentially with time. Moreover, the total 
change in azimuth angle is less than l80° in each case. Finally, the 
ferry mass loss is considerably closer to the : deal value in each of 
these cases than in the corresponding case for £ =0-7. 

For those cases in figure 6 for which £ - 0 . 7 , the distance, 
azimuth angle, speed, and acceleration all change much less regularly. 
Furthermore, the total change in azimuth angle is in each case much 


r— I Lf\ Lf\ 
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larger than for the corresponding ones for which £ = 1.0 and £ =1.3. 

In particular, for £ = 0.7 and v 2 = 0, collision must occur unless 
it is averted by applying more thrust than the control system requires. 
(See fig. 6, parts (a) and (b), where collision corresponds to an abrupt 
change of azimuth angle. ) 

In the interest of brevity, the variations of the direction angles 
of the vectors x and x have not been included in figure 6. They are 
readily found, if needed, from the appropriate equations in a manner 
similar to that used to compute the azimuth angle 0 as a function of 
time. 

Although figure 6 displays certain disadvantages of choosing £ = 0. 7 
it fails to give a clear comparison of the remaining cases, £=1.0 and 
£ = 1. 3> with one another. For this reason, the ratio of the ferry accel- 
eration to its initial value is plotted in figure 7 both as a function 
of the ratio of distance to initial distance and as a function of the 
ratio of speed to initial speed. Figure 7 shows that, in most cases of 
practical interest, the range of thrust acceleration required to reduce 
either distance or speed to a specified fraction of the initial value is 
substantially less for £ = 1.0 than for £ =1.3* 

The application of the results of figures 6 and 7 to a specific 
rendezvous problem necessitates more detailed considerations than the 
broad features which have been mentioned. When the basic rendezvous 
requirement is taken into account, together with the rate at which dis- 
tance and speed are reduced, the smoothness of change of the variables 
concerned, and the mass expended in propulsion, the choice of £ = 1.0 
seems, in the absence of extreme design requirements, to represent a 
reasonable starting point for examining the applicability of the present 
system. Thus, the one-dimensional analysis presented in the preceding 
section is particularly significant. 


EFFECT OF GRAVITATION 


The effect of gravitation will be examined for cases in which the 
primary attraction is that of the earth. Let | and if be the posi- 
tion vectors of satellite and ferry, respectively, in a rectangular 
Cartesian coordinate system with origin at the center of the earth and 
axes which correspond in direction and sense to the unit vectors of 
equation (l). The equation of motion of the satellite is, neglecting 
oblateness, and so forth, 

p— > — > 

d s GM i 

at 2 | 2 £ 


( 23 ) 
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where M denotes the mass of the earth and G denotes Newton's universal 
gravitational constant. Similarly, the ferry's equation of motion is 


d 2 rf = GM f F 

dt 2 " " H 2 t, m 


(24) 


Since 


T) = I + X 


(25) 


it follows that 


d 2 x _ d 2 ^ d 2 ^ 

dt 2 dt 2 dt 2 


( 26 ) 


Equations (23) to (26) yield the equation of relative motion 


* 2 | . I + gm/X _ A+JL 

-2 m I 1 3 i_* i5 


dt" 


4 + x 


(27) 


For the system of thrust application expressed in equation ( 3 ) > 
equation ( 27 ) "becomes 


dfx 

dt 2 


-2^ 0 


dx 

dt 


- O) 0 2 x + GMl^r 


— » — > 

I + x 


£ + X I 


( 28 ) 


In terms of its components, 


I = Sji + t 2 i + 


(29) 


— ) — ^ ) 1 k 

where the unit vectors i, j, k are the sane as those in equation (1), 


Numerical integrations of equation (28) lave been carried out for 
a circular satellite orbit described by the parametric equations 

, 2mt ^ 

£i = i sin — 


2nt v 
= -£ cos — \ 


(50) 




J 


vji oj vji h tr* 



where the orbital period T is given by 


JL 

2rt 



( 31 ) 


The radius selected was £ = 2.354 X 10^ feet, which corresponds to an 
altitude of approximately 2.64 X 10^ feet or 500 statute miles. Con- 
sideration was restricted to in-plane cases, that is, x^ = 0. Con- 
sequently, the ferry's position is given by the range (distance) 


R 



+ x, 




and the angle 


0 



(32) 


(33) 


The values 0.7 and 1.0 were used for the parameter £ and values 0.005> 
0.01, and 0.04 per second for the parameter odq. 

The results of the numerical computations are presented in fig- 
ures 8 to 12. The dashed curves, for which gravitation was neglected, 
are presented for comparison. The agreement is good for these sets of 
initial conditions. Thus, the simplified approach (gravitational effects 
being neglected) which was adopted earlier in this paper is justified 
for the approximate analysis of comparable cases. 


CONCLUDING REMARKS 


The linear system which has been considered for variable ferry- 
vehicle thrust control provides most of the relative-speed reduction 
early in the terminal phase. Position and velocity measurements rela- 
tive to a stable platform or its equivalent are required. 

The system 1 s fuel economy was found to be good for values of the 
control parameter £ of 1.0 and 1.3* In particular, it approaches 
the optimum as the normal component V *2 of the initial relative veloc- 
ity approaches zero. 

Thrust which is variable in magnitude as well as in direction is 
required in this system. Ways of meeting this requirement have been 
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indicated, but the degree of success in a given situation may depend 
upon the early development of rocket motors which can be used in this 
application. 

The required range of variability of thrust acceleration is sub- 
stantially less in most cases for £= 1.0 than for £ =1.3* Since 
a value of 0.7 for £ was found to be undesirable in other respects, 
the choice of a value of £ in the neighborhood of 1.0 is favorably 
indicated by the results which have been presented. The selection of 
the remaining control parameter o^q should be consistent with the design 
initial conditions . 


Langley Research Center, 

National Aeronautics and Space Administration, 
Langley Field, Va. , June 29, 19&1. 
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(a) Thrust program for a system of type (b). 



(b) Approximation by system of type (a) for which thrust interval is 
fixed, but number of rockets fired is variable. 



(c) Approximation by system of type (b) for which thrust magnitude is 

fixed, but duration is variable. 



Opposition angle 
(Variable) 


(d) Variable resultant thrust achieved by combination of constant- thrust 
motors with variable opposition angle. (Generally poor fuel economy. ) 

Figure 1.- Schematic illustration of methods for approaching a system of 

type (b) by systems of type (a). 
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Igure 2.- Position and velocity of ferry in a nonrotating (star- 
oriented) coordinate system with origin at satellite. 
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(a) Ferry launched from earth* 

Figure 3* “ Schematic illustration of desired situation at start of 
terminal- phase thrusting (for time t = 0). 








0 4 8 12 16 


Time, I 

(a) V x = -1.0; V 2 = 0. 

Figure 6.- Graphs of the dimensionless variables relative position (dis- 
tance and angle), speed, acceleration, and mass variable against time 
for various values of £ and initial relative velocity Vn, 1 ^- 












Rotio of ferry-to-satellite 
distance to initial value 


Ratio of ferry s relative 
speed to initial value 


( a ) Vl = -0.5. 

Figure 7-- Graphs of thrust acceleration against distance and speed, all 
referred to their initial (t = 0) values, for £ = 1.0 and £ =1.5* 
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sec Time, sec 

(e) x x = 100,000 feet; (f) x = 100,000 feet; 

i]_ = -500 feet per second. x 1 - -1,000 feet per second. 




Time, sec 


(g) x-l = 1,000,000 feet; 

= -1,000 feet per second. 



Time, sec 


(h) x L = 1,000,000 feet; 

= -2,000 feet per second. 


Figure 8.- Concluded. 
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Time, sec 

(a) x-^ = 1,000 feet; 

x^ = -20 feet per second. 




Time, sec 

(b) x-^ = 1,000 feet; 

x^ = -100 feet per second. 




Time, sec Time, sec 


(c) x^ = 10,000 feet; (d) x^ = 10,000 feet; 

= -100 feet per second. x^ = -500 feet per second. 

Figure 9- - The effect of gravitation on the ferry's motion relative to 
satellite in a circular orbit at an altitude of 500 statute miles for 
£ = 0.7 and gdq = 0.0*1 per second. 
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0 800 1600 2400 3200 

Time, sec 


(e) x-^ = 100,000 feet; 

= -500 feet per second. 



Time, sec 


(g) x 1 = 1,000,000 feet; 

x^ = -1,000 feet per second. 


Figure 10. - 



0 800 1600 2400 3200 

Time, sec 


(f) = 100,000 feet; 

x^ = -1,000 feet per second. 



(h) x-^ = 1,000,000 feet; 

x-j^ = -2,000 feet per second. 


i. 
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Time, sec Time, sec 

(a) x-^ = 1,000 feet; (b) x-^ = 1,000 feet; 

x-^ = -20 feet per second. x-^ = -100 feet per second. 



Time, sec Time, sec 


(c) x x = 10,000 feet; (d) xp = 10,000 feet; 

Xp = -100 feet per second. Xp = -500 feet per second. 

Figure 11.- The effect of gravitation on the ferry's motion relative to 
a satellite in a circular orbit at an altitude of 500 statute miles 
for £ =1.0 and = 0.04 per second. 



k2 



(e) x-j_ = 100,000 feet; 

x-^ = -50O feet per second. 



(f) xj = 100,000 feet; 

X]_ = -1,000 feet per second. 




Time, sec Time, sec 


(g) x 1 = 1,000,000 feet; (h) x x = 1,000,000 feet; 

x = -1,000 feet per second. x ± = -2,000 feet per second. 


Figure 11. - Conclude 1. 
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0 400 800 1200 1600 

Time, sec 


(a) x 1 = 100,000 feet; 

xj_ = -500 feet per second; 

x 2 = 

X 2 = -250 feet per second; 

i = 0 . 7 . 



0 400 800 1200 1600 

Time, sec 


(c) x-j_ - 100,000 feet; 

= -500 feet per second; 
x 2 = 0; 

x 2 = -250 feet per second; 

g = 1.0. 



O 400 800 1200 1600 


T i me, sec 

(b) Xj = 100,000 feet; 

x^ = -500 feet per second; 
x 2 = 50,000 feet; x 2 = 0; 

<L = 0 . 7 . 



O 400 800 1200 1600 


Time, sec 

(d) = 100,000 feet; 

x-^ = -500 feet per second; 
x^ = 50,000 feet; 

X 2 = 0; £ = 1. 0. 


Figure 12.- The effect of gravitation on the ferry’s motion relative to 
a satellite in a circular orbit at an altitude of 500 statute miles 
for gdq = 0. 01 per second. 
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